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Abstract. MASSCLEAN is a new, sophisticated and robust stellar cluster image and photom- 
etry simulation package. This package is able to create color-magnitude diagrams and standard 
FITS images in any of the traditional optical and near-infrared bands based on cluster character- 
istics input by the user, including but not limited to distance, age, mass, radius and extinction. 
At the limit of very distant, unresolved clusters, we have checked the integrated colors created in 
MASSCLEAN against those from other simple stellar population (SSP) models with consistent 
results. Because the algorithm populates the cluster with a discrete number of tenable stars, it 
can be used as part of a Monte Carlo Method to derive the probabilistic range of characteris- 
tics (integrated colors, for example) consistent with a given cluster mass and age. We present 
the first ever mass dependent integrated colors as a function of age, derived from over 100,000 
Monte Carlo runs, which can be used to improve the current age determination methods for 
stellar clusters. 
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1. Introduction 

Star clusters are the building blocks of galaxies. In order to analyze the known pop- 
ulation of star clusters, to search for new clusters, and to derive the selection effects of 
current NIR survey searches for Milky Way clusters we created a simulation package. 

The MASSCLEAN packag^J (Popescu & Hanson 2009; Hanson & Popescu 2007, 2008, 
2009) is built using numerous well established theoretical and empirical models for stars 
and stellar clusters, starting with the Kroupa-Salpeter IMF for stellar mass distribution 
(Kroupa 2002; Salpeter 1955). The user has the option to choose between two stellar 
evolutionary models: the Geneva Models (Lejeune & Schaerer 2001) or the Padova Models 
(Marigo et al. 2008). The extinction model is based on Cardelli, Clayton and Mathis 
(1989), and the King Model is used for the spatial distribution (King 1962). An optional 
stellar field can be generated using the "SKY" Model (Bahcal & Soneira 1984). MASSCLEAN 
is built in a modular way and is able to combine observational data and simulations using 
an unprecedented level of accuracy. 

2. Images — Galactic and Extragalactic Simulations 

MASSCLEAN generates FITS images, using SkyMaker (Bertin 2009). The simulated J- 
band image for NGC 3603 is presented in Fig.[T](b), and the 2MASyil image is presented 
in Fig.n(a). 

t The package is freely available under GNU General Public License at: 
http : //www . physics . uc . edu/~bogdan/mass clean/ 

^ This publication makes use of data products from the Two Micron All Sky Survey, which is a joint 
project of the University of Massachusetts and the Infrared Processing and Analysis Center/California Insti- 
tute of Technology, funded by the National Aeronautics and Space Administration and the National Science 
Foundation. 
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Real Image (2MASS) (b) Simulated Image 

Figure 1. ATGC 3603 J-Band Images (4.4 x 4.4 arcmin). 

In the Fig. [2] we present a simulated 10^ Mq cluster at the distance of M31 and the 
resolution of HST ACS. The images show the clusters at ages of 1 and 100 million years. 




(a) log{age/yr) = 6.0 (b) log{age/yr) = 8.0 

Figure 2. 10^ Mq simulated cluster at the distance of M31 and the resolution of HST ACS. 

We are using MASSCLEAN to help us derive the selection effects of current NIR survey 
searches for Milky Way clusters (Hanson et al. 2009; Hanson et al., in prep.). As a proof 
of concept, we present an example image simulation of 45 randomly generated clusters 
of various age, mass, concentration, and extinction. The F-band image is shown in Fig. [3] 
(a), and the i^T-band image is shown in Fig.[3](b) (note, no stellar field is included). The 
positions of the clusters along with the scaled angular size of the clusters' core is shown 
in Fig.[3] (c). 
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(a) V-Band Image (b) K-Band Image 

Figure 3. 45 simulated clusters 



(c) Positions 



3. Integrated Colors and Ages of Clusters 

MASSCLEAN follows the evolution of every star in the cluster, providing full pho- 
tometry at all ages. We have tested that at the infinite mass limit the integrated colors 
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are consistent with standard SSP models. A log {age /yr) step of 0.05 has been used to 
allow comparison with all other models. The integrated colors generated by MASSCLEAN 
using Padova 2008 evolutionary models for solar metalicity are presented in Fig.[H (a). 
The results from other SSP codes (Marigo et al. 2008; Bruzual & Chariot 2003; Maraston 
2005; Kotulla et al. 2009) are also displayed, along with the Milky Way clusters from 
Hancock et al. 2008. The MASSCLEAN integrated colors follow the Marigo et al. 2008 since 
they both use the same evolutionary models for their input. 
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(a) Comparison with modern SSP codes (b) Mass dependent integrated colors 
Figure 4. MASSCLEAN integrated colors 



Of course, real star clusters have a finite mass. They do not have fully populated 
isochrones at all ages and their integrated colors may not always correspond to the 
ones computed in the infinite mass limit (Bruzual 2002; Cerviho & Luridiana 2004; 
Cerviho & Luridiana 2006; Fagiolini et al. 2007; Lancon & Mouhcine 2000; Lancon & 
Fouesneau 2009). The dispersion of integrated colors (Icr) due to different values of mass 
and fiuctuations in the IMF is presented in Fig. [4] (b). {B — V)o vs. log{age/yr) and 
(U — B)o vs. log {age /yr) are presented for a few values of cluster's mass. The step in 
log {age /yr) is 0.01. The color dispersion for clusters with mass in the 2 x 10^ Mq - 
2.5 X 10^ Mq range is based on about 100,000 Monte Carlo runs (which corresponds to 
about 15 million cluster models). We present for the first time mass dependent integrated 
colors for clusters. The entire set of mass dependent integrated colors using both Padova 
and Geneva models, and different metalicities will be publicly available soon (Popescu & 
Hanson, in prep.). 

Age determination for clusters using integrated colors is usually based only on the 
infinite mass limit. Using our mass dependent integrated colors and magnitudes, we 
expect to provide a more realistic model. As an initial demonstration, we present new 
age determination based on MASSCLEAN integrated colors for numerous Milky Way clusters 
in Fig.O (Popescu & Hanson, in prep). 
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(a) U-B colors as a function of age. (b) B-V colors as a function of age. 



Figure 5. Real clusters, with masses of just a few thousand solar mass, are expected to be 
redder than SSP models predict using an 'infinite mass assumption' particularly in the age 
range 10^ — 10^. We are testing this result using integrated colors of Milky Way clusters with 
well constrained ages (Popescu & Hanson 2009, in prep.). 

4. Summary 

We have developed a sophisticated cluster image and photometry simulation package 
which we are using to test against current search algorithms. We predict that clusters 
over certain age ranges (particularly when the red supergiant phase is just starting out) 
and clusters with low stellar densities (even though possibly of very high total mass) are 
being missed by the current searches, particularly when combined with high extinction. 
Our mass dependent integrated colors for clusters should provide a better description for 
known clusters and a more accurate age determination. 
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